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1.  Introduction 


™The  analysis  of  steady  plane  deflagration  vaves  invariably  starts  vith  the 
combustion  approximation  where  it  is  assumed  that  the  Mach  number,  i.e.  the  flame 
speed  divided  by  a  characteristic  sound  speed,  is  vanishingly  small.  The  momentum 
equation  then  implies  that  the  pressure  is  nearly  constant  while  the  thermal 
and  mechanical  descriptions  of  the  wave  decouple,  so  that  the  task  of  solving 
for  the  structure  is  greatly  simplified.  Even  then  explicit  formulas  can  only 
be  obtained  in  the  limit  of  large  activation  energy. 

We  are  currently  interested  in  describing,  by  means  of  activation-energy 
asymptotics  as  far  as  possible,  the  transition  from  deflagration  to  detonation 
in  gases. *  One  of  the  first  steps  in  such  a  theory  is  to  analyze  deflagration 
vaves  whose  Mach  numbers  are  not  vanishingly  small.  Pressure  variations  cannot 
be  neglected  and  hence  the  momentum  equation  must  be  retained  in  the  descrip¬ 
tion  of  the  structure.  We  will  show,  in  the  present  paper,  that  the  method  of 
activation-energy  asymptotics  gives  an  analytic  description  of  these  fast 
deflagrations,  i.e.  deflagrations  travelling  at  Bpeeds  greater  than  those  Justi¬ 
fying  the  use  of  the  combustion  approximation.  In  addition  we  examine  the  limit 
of  vanishingly  small  Mach  number  to  shed  light  on  the  nature  of  the  combustion 
approximation, 

The  main  result  is  that  the  Mach  number  of  a  fast  deflagration  wave  is 
determined  (implicitly)  by  the  flame  temperature ,  a  parameter  of  the  mixture. 

But  an  interesting  feature  is  a  secondary  reaction  over  an  exponentially  long 
tail  in  which  the  small  amount  of  reactant  escaping  the  initial  flame  is  finally 
consumed.  In  addition,  for  small  Mach  numbers  we  find  a  hydrodynamic  adjustment 
layer  behind  the  flame  and  another  possible  wave^neither  of  which  can  be  detected  by 
the  combustion  approximation^  In  short,  we  have  extended  the  classical  work 
of  Bush  and  Fendell  (1970)  to  finite  Mach  numbers,  and  uncovered  the  existence 
of  very  slow  deflagration  vaves. 
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2.  The  governing  equations 

The  governing  equations  of  a  plane,  steady  deflagration  express  the  balances 
of  mass,  momentum,  energy  and  species.  In  deriving  the  equations  used  here  cer¬ 
tain  assumptions  are  made,  the  most  important  being:  steadiness, one-step  reac¬ 
tion,  Arrhenius  kinetics,  ideal  gas,  Newtonian  fluid.  Pick's  diffusion  lav, equal 
specific  heats  and  constant  material  properties.  See  Buckmaster  &  Ludford  (l98l). 

In  what  follows  p,v,T  and  Y  are  respectively  the  dimensionless  density, 
gas  velocity,  temperature  and  mass  fraction  of  the  deficient  reactant  (if  there 
is  more  than  one).  The  density  and  temperature  in  the  quiescent  conditions 
upstream  are  taken  as  units  so  that  the  dependent  variables  p,v, T,Y  tend  to 

l,0,l,Y_w  as  s  -*•  -«,  where  s  denotes  distance  from  the  flame.  The 
distance  unit  is  the  preheat  thickness  X/c^M,  where  x  and  c^  are  respectively 
the  thermal  conductivity  and  specific  heat  of  the  fluid;  and  M  is  the  mass  flux, 
the  velocity  unit  being  the  flame  speed.  The  other  important  parameters  that 
appear  are  defined  in  Table  1. 

The  Prandlt  and  Levis  numbers  will  both  be  set  equal  to  one,  solely  to 
simplify  our  discussion.  Prandlt  and  Lewis  numbers  different  from  unity  8LT€ 
discussed  in  Appendix  I. 

The  equations  of  steady  plane  combustion  can  then  be  shown  to  reduce  to 


Y2^dV/ds  -  YI^(V-l)  +  (TV_1-l),  (l) 

o  o 

dt/ds  -  d2t/ds2  +  oAYe“0/T,  (2) 

T  -  t  -  (Y-l)J^V2 /2,  (3) 

o 

t  +  aY  »  1  +  8+  (Y-l)f^/2  with  8  ■  OY  ,  (4) 

o  — • 

pV  ■  1  -with  V  ■  v  ♦  1.  (5) 
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The  perfect-gas  law  serves  to  define  the  pressure,  which  has  been  eliminated. 

In  the  present  formulation,  V  and  x  will  he  the  basic  variables,  the 
temperature  T  and  mass  fraction  Y  serving  as  auxiliary  variables  defined  by 
(3)  and  (4)  respectively.  Thus  equations  (l)  and  (2)  are  to  be  solved  under 
the  conditions  that  V,x  -*■  l,l+(  Y-l)i^/2  as  s  -•  and  that  the  solution 
is  bounded  as  s  -*•  +».  Equation  (3)  Just  defines  the  density.  The  singular 
nature  of  equation  (l)  as  Mq  -*  0  leads  to  the  adjustment  layer  mentioned 
in  the  Introduction  (cf.  Secs.  5  and  6). 

3.  Activation-energy  asymptotics 

We  seek  a  solution  of  equations  (l)-(U)  in  the  distinguished  limit 

D  »  C02exp(6/T#)  with  6  ■*■»  .  (6) 

Here  the  constants  C  and  T„  are  supposed  given;  their  determination 
for  an  actual  mixture  is  considered  in  Sec.  8.  The  connection  between  the 
wave  speed  and  the  mass  flux  H  is  made  explicit  by  writing 

M  *  mM  (7) 

o 

(m  is  the  so-called  acoustic  impedance).  Thus  the  reaction  term  in 
equation  (2)  appears  as 

AYe"6/>T  «  C02Yexp[6(l/T#-l/T>]/m2^  (8) 

As  0  •,  there  are  three  possible  regimes.  In  regions  where  T  <  T„, 

the  reaction  term  is  exponentially  small  and  hence  negligible  compared  to 
algebraic  perturbations;  the  chemistry  is  frozen.  In  regions  where  T  >  T#, 


the  reaction  term  is  unbounded  unless  Y  is  zero  to  all  order  in  0”^ ; 
equilibrium  prevails.  Finally,  the  regions  in  vhich  T  -  T#  »  oCd’^)  are 
flame  sheets,  found  to  be  characterized  by  a  reaction/diffusion  balance. 
Equilibrium  regions  will  in  fact  only  occur  in  the  limit  of  the  combustion 
approximation. 


4.  Fast  deflagration  waves 

We  now  determine  the  structure  of  fast  deflagration  waves.  Bie  flame  sheet 
is  located  at  s  *  0  and  on  either  side  the  chemistry  is  frozen,  as  we  shall 
see.  Expanding  the  variables  in  the  form 


u  *  uq  +  e"1^  +•••  (9) 

and  assuming  Y  is  zero  to  leading  order  in  the  burnt  region  shows  that 

t  -  1  +  0eS  +  (y-1)J^/2,  Y  «  Y  (l-e8)  for  s  <  0  ,  (10) 

O  O  O  -m 

and 

t  -  1  +  8  +  (Y-l)J^/2,  Yrt  -  0  for  s  >  0  .  til) 


The  solution  for  Vq  requires  some  discussion.  On  both  sides  of  the 

flame  sheet  VQ  is  governed  by  (l)  where  V  and  T  are  taken  to  leading 

order  and  T  is  defined  in  terms  of  V  and  t  by  (3).  For  s  >  0  the 
o  00 

fixed  points  of  equation  (l)  are 

V+  -  [(Y^+l)  ±  /(l-M?)2  -  2(Y+l)tfW(Y+l)tf;  Cl2) 

z  O  O  OO 

which  shows  that  the  restriction 

(1-I^)2/2(y+1)^  >  8  (13) 

must  be  placed  on  Mq  for  a  given  heat  release  8  .  Equality  corresponds 
to  the  Chapman- Jouget  deflagration,  so  that  the  present  theory  will  yield  Mach 


numbers  in  the  range 
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0  < 


M2 

o 


[l  ♦  (r*D  B-  /U+(t*i)b)2-i]. 


ilk) 


Integrating  (l)  for 


s  >  0  shows  that 


2y  1 

Y+l  (V+-V_) 


in 


s  +  c 


(15) 


where  c  is  a  constant.  Since  V+  -  V  is  positive,  VQ  V_  as  s  -*■  +«#, 
which  is  appropriate  for  weak  deflagrations. 

For  s  <0,  we  cannot  integrate  equation  (l)  analytically  without  further 

g 

assumption  because  Tq  contains  e  .  Nevertheless  its  integration  with  V  ®  1 

for  s  a  -«  determines  the  value  V  ■=  VM  (M2)  <  V  at  s  *»  0.  Determination 

o  o  - 

follows  from  the  saddle-point  nature  of  the  cold- boundary  point  in  the  x,V-plane; 

the  inequality  results  from  dV  /ds2  O  for  8  >  0,  not  a  simple  result  to  verify. 

o 

As  a  consequence  Tq  decreases  from  its  value  T„  at  the  flame  because  tq  is 
constant . ) 

Since  V„  is  determined  by  the  integration  for  s  <0,  the  constant  c 
in  (15)  is  fixed  and 


*  tq{o)  «  1  +  a  +  (y-i)m^(i-vJ)/2. 


(16) 


I’hus  the  flame  temperature  is  determined  as  a  function  of  the 

flame  speed,  i.e.  If  T#  is  a  monotonically  decreasing  function  of  1^, 

which  is  usually  found  to  be  the  case  for  reasonable  values  of  a»Y»Pr*  "then 

there  is  a  unique  M2  for  the  given  T#  in  the  definition  (6)  of  D.  A  typical 

plot  of  vs.  T#  is  shown  in  Fig.  la. 

Note  that  no  discussion  of  higher-order  terms  or  structure  of  the  flame 

sheet  is  necessary  to  determine  the  flame  speed.  Only  the  minimal  assumption 

of  continuity  across  the  flame  sheet  is  made.  The  analysis  is  distinctly  different 

from  that  for  M  small  (using  the  combustion  approximation),  where  a  discussion 
0 

of  the  structure  (and  hence  higher-order  terms)  is  needed  to  determine  Mq. 
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On  either  side  of  the  flame  sheet  the  perturbations  of  t  and  Y  are 

t1  »T_es  for  s  <  0,  *  T+  for  8  K  0  8114  \  *  -Tj/a»  (17) 

while  satisfies 

YI^dV,/ds  =  V.^-T  /V2)  +  t/Vft  for  all  s  4  0.  (l8) 

ol  looo  10 

The  integration  constants  t+  are  determined  by  the  flame-sheet  structure,  as 
we  shall  see  next.  Note  that  t+  is  not  zero:  an  0(e""^)  amount  of  reactant 
escapes  burning  and  passes  downstream  where,  over  exponentially  large  distances, 
it  is  finally  consumed.  The  process  is  described  in  Appendix  II. 

To  investigate  the  structure,  let 

z  ■  8s  (19) 


and  consider  the  expansions 

T-  1+8+(y-1)»^/2  ♦  e“1T1(z)+-*- ,v  -  V#  +  e“1V1(t)+... ,  Y  *  6~1Y1(z)+...  (20) 

T  ■  T#  +  0-1T1(z)+. . . 

Then  we  easily  find 

»  «z  +  v  (0)  with  6-  V#  -  1  +  (WDyM2,  (21) 

where  V^(0)  is  at  this  stage  an  unidentified  constant.  However,  it  follows 
by  matching  that  in  the  outer  expansions  is  continuous  across  the  flame 

sheet,  so  that  V^(0)  is  in  fact  the  common  value  there  (hence  the  notation). 

The  solution  of  equation  (l8)  is  determined,  once  t+  are  found,  by  the  upstream 
condition  V^(— «•)  ■  0.  Clearly  the  upstream  solution,  and  hence  the  constant 
V^(o)  in  (21),  depends  on  t  while  the  downstream  solution  depends  on  both 
t+  and  t_. 
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The  -problem  now  becomes 
T  /T2 

d2f,/dz2  =  AT.e  1  *  with  T,  ■  t,  -  ( Y-l [  4z+V,  (0>) 
1  1  11  o  w  l 


and  A 


C/m2!^, 


■  gz  +  T  +  o(l)  as  Z  ■+  -« ,  *  T+  +  o(l)  as  z  ■+•  +oo 


where  the  boundary  conditions  are  obtained  by  matching.  A  transformation 
reduces  (22)  to  a  structure  problem  discussed  by  Linan  [197*0 .  Let 


o  =  -tj/T#,  w  *  (y-1)M^V#6/s, 

C  “  ez/T 2  +  6V1(0)/6T2  +  (l/w)jm(B2/2AT*}; 


2  2 

2d  a/d?  »  oexp{-o-w^}. 


o  *  -C-c_/T^  +  V^Oj/gT2  +  (l/w)An{g2/2AT^}  +  o(l)  as  c+-«A(2*0 

o  -  -t+/T2  +  o(1)  as  +oo.  I 

Lilian  found  (numerically)  that,  whenever  w  is  less  than  1/2,  a 
solution  of  the  equation  (2*t&)  is  uniquely  determined  by  the  weaker  boundary 


conditions 


do/d£  -+-1  as  5  ■+•-08  ,  do/d^  ■+■  0  as  5  +o+. 


By  adjusting  the  (finite)  starting  point  for  the  left  boundary  condition 
he  was  able  to  satisfy  the  right  boundary  condition,  and  thereby  approximate 
the  parameters 


o_(v)  ®  lim  [a* c),  0  (w)  ■  lim  o* 

5-+— 00  £-»+oo 


(26) 
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For  0  <  w  <  1/2,  as  is  the  case  here,  o+  was  found  to  be  positive  while 
o  -  1.3UU,  o+  -*•  0  as  v  0,  [Proof  that  w  is  less  than  1/2  has  been  given 
by  Lu  4  Ludford  (l98l).]  It  follows  that  for  our  problem 

-t_/T2  +  V1(0)B/6T^  +  W'snU^AT^}  »  a_,  -t+/T*  «  o+  (27) 


are  therf^ determined. 

To  complete  the  determination  of  r+  we  must  calculate  V^to)  in 
terms  of  t  ,  which  requires  the  equation  (l8)  for  in  s  <  0  to  be 

integrated  tinder  the  condition  V^(-“)  -  0.  The  solution  is 

V  »  T  fl(s),  with  fl(s)  ®  f  exp  ds'  (2E 

Y VT  -»  YiTV  (s') 


where 


2 

f(s)  =  -V  f  (rfVT  /V*)ds': 
YM2  o  000 


and  hence 


v]L(o)  -  T_fi(0). 


We  now  have  asymptotic  approximations  to  order  0  in  the  three 
regions,  for  any  set  of  parameter  values  C,T# ,m, |5,y,Y  To  summarize^ 
these  approximations  are  determined  as  follows.  First  t  ,  and  hence 
Yq  j  is  determined  by  solving  the  reactionless  equation  (2)  on  either 
side  of  the  flame  sheet  assuming  continuity  across  the  sheet  and  Yq  *  0 
for  s  >  0.  The  solution  is  given  by  equations  (10)  and  (ll).  Vq  is 
then  determined  by  integrating  equations  (l)  to  leading  order,  numerically 
for  s  <  0  and  analytically  [equation  (15)1  for  s  >  0.  The  temperature 


follows  from 


^uatio|  (3), 


depending  in  particular  on  Mq.  Setting 


Tq(0)  *  T#  then  fixes  the  flame  speed  (represented  by  M^)  as  a  function 
of  the  -PIclvy\£  temperature  T#,  according  to  equation  (l6).  Pig,  2a  is 
the  result  of  carrying  out  these  steps  for  one  particular  set  of  parameters. 

The  approximation  to  the  next  higher  order  is  obtained  by  solving 
for  and  in  the  outer  regions  (s  £  0),  S®6.  equations  (17) 

and  (28),  and  determining  the  corresponding  constants  x+,V^(0)  through 
the  relations  (27)  and  (29).  These  relations  are  obtained  from  the  structure 
of  the  flame  sheet  and  explicit  integration  of  the  equation  for  in 

s  <  0,  respectively. 

Explicit  analytical  results  can  be  obtained  when  the  heat  release 
is  small,  in  which  case  all  variables  stay  close  to  their  values  upstream. 
Note  that  this  limit  can  be  taken  for  any  value  of  Mq,  i.e.  it  is  quite 
independent  of  the  combustion  approximation  Mq  •+•  0. 


Small  heat  release:  6  «  1 

The  smallness  of  the  heat  release  may  be  due  either  to  a  small  heat  of 
reaction  or  to  a  small  amount  of  reactant.  It  enables  analytical  expressions 
to  be  developed  for  the  results  in  Sec.  4  that  otherwise  have  to  be  derived 
by  numerical  integration. 

We  write 

Vq  *  1  +  gV'  +•••  ,  tc  »  1  +  (y-l)^/2  +gT'  +•••  (30) 

so  that 

2 

Y o  =  &(l-x' w-IT,  -  1  +  B[t'-(Y-1)MoV']  +...  .  (31) 
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YM^V'/ds  -  -(l-M^V*  +  t'. 


The  solution  which  vanishes  at  s  =  -*>  and  is  continuous  at  s  *  0  is 


f 

[l+(Y-l)j^]""^eS  for  s  <  0, 
o 

[l-M?]“1{l-Y^tl+(Y-l)^r1e  ^  for  s  >  0  . 


Note  that  when  Mq  is  small  a  rapid  hydrodynamic  adjustment  takes  place 
behind  the  flame  sheet,  given  by  the  exponential  term  in  (34  b).  The  combustion 
approximation  takes  the  thickness  of  this  Mach  layer  [which  is  a  consequence 
of  the  singular  nature  of  the  differential  equation  (33)  as  Mq  -*•  0]  to  be  zero. 

Consider  the  distinguished  limit 


T#  *  1  +  Bt#,  C  *  B  C  with  g  0. 


For  t#  in  an  appropriate  interval,  T#  is  attainable,  i.e.  it  lies 
between  the  upstream  temperature  and  the  adiabatic  flame  temperature;  and 
t  ,  obtained  from  equation  (27  a),  remains  bounded  as  g  0.  As 
explained  in  Sec.  4,  in  order  to  determine  the  wave  speed,  we  calculate 
the  temperature  at  s  *  0  and  set  it  equal  to  T#.  Here  the  relation 
between  T#  and  becomes  explicit ,  namely 


t#  =  1/[1+(y-1)^3 


In  particular,  as  T#  approaches  the  adiabatic  flame  temperature  Tft,  i.e. 

t«  -*•  1,  the  flame  speed  vanishes,  i.e.  M2  -►  0.  The  restriction  (13)  on 
w  o 

the  heat  release  is  not  violated  unless  becomes  close  to  M  ■  1  +  Otg1^2), 

o  OCT 

thus  t#  varies  between  1  and  1/Y  <  1. 
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To  complete  the  analysis  ve  calculate  the  perturbations  explicitly. 

Details  are  omitted;  suffice  it  to  say  that,  to  leading  order  in  $,  we  find 

d  *  0/Y,w  =  [(Y-1)/Y]l£,  fiC s )  «  [l/(l+(  Y-l)M^)]eB  (37) 


so  that 

V_(0)  =  t  /[1+(Y-1)J^], 
l  -  o 


(38) 


where,  from  equation  (27  a), 


t 


l+tY-ljM2 

o 

( Y-i) (l-M^ ) 


(39) 


In  short,  the  solution  to  0( 8-1) ,  both  Inside  and  outside  the  flame  sheet, 
is  explicitly  determined  once  C  and  t#  are  specified.  The  flame  speed  is 
determined  by  t#  alone,  while  2  only  enters  at  the  perturbation.  This 
means  that  changing  C  affects  the  structure  of  the  flame (corresponding  to 
a  slight  translation^and  hence  the  perturbation  it  produces  upstream  and 
downstream,  but  not  its  velocity. 


6.  The  Combustion  Approximation 

The  limit  M  -*■  0  must  be  examined  carefully.  But,  before  doing  so,  we 
o 

shall  review  the  basic  results  of  activation-energy  asymptotics  in  the 
combustion  approximation. 

The  combustion  approximation  is  the  limit  Mq  0  vith  0  fixed.  The 
activation-energy  limit  Q  •*  <»  is  then  used  to  analyze  the  resulting 

equations.  Thus  lim  lim  is  obtained  in  contrast  to  lim  lim  which  will 

0-h»  M  -*0  M  -*0  0-*» 

o  o 

be  obtained  from  the  above  analysis.  The  pre-exponential  factor  C  in  the 
Damkohler  number  is  assumed  to  be  such  that  A  ■  C/m^M?  is  0(l). 


i 

i 

1 

1 
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To  all  orders  in  equations  (1-5)  with  M  *  0  are 

o 


correct 


to  leading  order  in  a  ^-expansion.  We  find  immediately 


T«x»V,  T+aY*l+e-T 


(1*0) 


to  all  orders  in  e""^;  and  the  leading  order  result 


To  -  < 


1  + 

Be8 

for 

8  <  0 

1  + 

S 

for 

8  >  0 

(1*1) 


Standard  flame-sheet  analysis  yields  the  structure  problem 

(22)  which,  with  t.,  ,  T#  replaced  by  T- ,T  and  M  ■  0,  becomes 

x  la  o 


9  /t 

d2^/dz2  -A^e  1  a  vith  A  *  </mV 

■  0z  +  r  +  o(l)  as  z  -*•  -oo,  *  x+  *  o(l)  as  z  +«. 


(1*2) 


Clearly,  t+  is  zero;  otherwise  the  differential  equation  is  not  consistent 
with  the  boundary  condition  at  z  ■  +».  A  first  integral 


i  T  l<JT 

(dT./dz)2  -  2T„A[l+(T,/T2-l)e  1  a] 
i  a  i  a 


(1*3) 


is  obtained  by  using  the  requirement  t+  ■  0.  The  boundary  condition 
at  z  ■  -oo  will  then  Ta e  satisfied  only  if 

p  li  p  P  9  li 

A  -  0  / 2r,  i.e.  DVTc  -  -  B  0  e  a/2TA. 


(1*1*) 


A  further  integration  using  the  same  condition  yields 

P2 

a 


n/T; 


8* 


-  t_  *  T2  j  {[l-(t-l)et)"1/2-l}dt  +  iv 


(1*5) 


\ 


Thus,  while  the  constant  t+  is  zero,  is  undetermined;  its  definition 

corresponds  to  locating  an  origin  for  the  flame  structure  to  0(0**^").  For 
the  structure  to  exist,  A  must  have  the  value  (44),  corresponding  to  a 
definite  distinguished  limit  for  D  as  0. 

Returning  to  the  fast  deflagration  with  0  ■*•  •  ,  we  now  let  +  0. 

First  note  from  the  discussion  in  See.  4  that  M  -*■  0  as  T«  -*•  T  :  the  limit 
of  vanishingly  small  flames  speeds  corresponds  to  fllme  temperatures  T# 
indefinitely  close  to  the  adiabatic  flame  temperature  T  . 

The  only  point  at  which  there  is  difficulty  in  the  limit  is  the  unbounded 

term 

w“1tn{02/2ATt}  (46) 

in  the  formula  (27).  If,  however,  we  set 

A  =  e2/2Tt,  i.e.  C  -  02ml£/2Ti,  (47) 

corresponding  to 

D  ■  B2ra2M^(T#)02exp(0/T#)/2Tt,  (48) 

the  term  disappears  and  no  difficulty  remains.  We  have  therefore  constructed 
a  family  of  asymptotic  solutions,  parametrized  by  T#,  which  is  uniformly 
valid  in  the  interval 

Vj<  t*I  V  (49) 

Moreover,  the  member  T#  *  T  of  this  family  is  one  of  the  solutions  obtained 
by  the  combustion  approximation. 
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Other  uniformly  valid  families  can  be  obtained  by  taking 
A  *  (62/2Ti)[l+M^(T#)f(T#)]  with  f  arbitrary.  This  would  add  a  factor 
1  ♦  ^(T,)f(T#)  to  the  right-hand  side  of  (48)  and  change  t_  by  an  amount 
depending  on  T#.  Because  of  the  exponential  in  the  formula  (48),  any  two  such 
families  produce  the  same  finite  D  for  values  of  T#  within  0(e  1)  of 
each  other.  The  corresponding  members  are  therefore  identical  except  for  T_, 
a  change  which  corresponds  to  translation  of  the  wave  by  an  0( e”1)  amount. 

The  two  members  corresponding  to  T#  ■  T  differ  in  the  same  manner, 

ft 

reflecting  the  Indeterminancy  noted  above  for  solutions  under  the  combustion 
approximat ion . 

When  Ta  is  excluded  from  the  T#-interval  a  much  larger  class  of  uniformly 
valid  families  of  solutions,  all  identical  under  translations,  can  be  obtained 
by  taking  A  to  be  an  arbitrary  function  of  T„,  again  because  of  the  exponential 
in  D.  Any  such  family  which  is  not  of  the  earlier  type  breaks  down  as 
T-  -*•  T  because  the  term  (46)  becomes  unbounded  in  the  limit,  suggesting 

W  ft 

that  there  are  other  solutions  when  M  is  small,  characterized  by 

o 

perturbations  that  are  large  compared  to  0  \ 

We  turn  now  to  these  solutions  of  our  general  equations,  finally 
demonstrating  that  they  do  not  survive  the  combustion  approximation.  It 
is  therefore  not  surprising  that  they  have  escaped  detection  before. 

7.  Very  Slow  Deflagrations 

We  now  consider  the  distinguished  limit 

D  *  C02exp(e/T#) ,  T#  ■  T&  -  g(e)t„  with  0  -►  «,  (50) 

where  g(o)  is  a  positive  guage  function  large  compared  to  0”^  and  t# 
is  0(l).  We  anticipate  M2  *  kg(o),  with  k  ■  0(1)  positive,  and  look 


for  expansions  of  the  form 


u  *  uQ  +  g(0)ug  +  ®”lul  +‘“  (51) 

upstream  and  downstream  of  the  flame.  Thus  the  leading-order  terms  are  given 
hy  formulas  (10)  and  (ll)  with  «  0  and  V  ■  T  *  t  •  The  g-perturbations 

O  0  0  0 

are  then  found  hy  assuming  continuity  of  T  and  t  (hut  not  V  )  from 

g  g  g 

one  side  to  the  other,  assumptions  that  can  he  justified  hy  considering  the 
structure  of  the  sheet.  Thus 


T 

g 


T 

g 


(Y-l)k/2  +  [Tgio-<Y-l)k/2]e8,  Vg  ■  Tg  -  (y-l)k//2  for  s  <  0, 

t  ,  V  «*  t  -  (y-l)kTf/2  +  kygT  for  s  >  0, 
g»  g  a  a 


(52 


where  r ^  is  an  as  yet  undetermined  constant.  As  expected  from  the  discussion 
in  Sec.  5,  V  is  discontinuous  across  the  flame  sheet  so  that  a  Mach  layer  is 

s 

required.  We  easily  find  the  modification 

V  -  T  -  (Y-l)kT*/2  +  YtaT  [l-exp(-s/ykT  g)]  for  s  >  0,  (53) 

g  9°  a  ”  a  a 

which  makes  the  necessary  adjustment  on  a  scale  s  ■  0(g),  much  larger  than 
that  of  the  flame  sheet.  The  0 ”^-perturhations  are  found  to  he 

*  t  e8  for  s<  0;  ®  T  +  for  s  >  0,  (5^) 

where  t±  have  to  he  determined. 

To  fix  the  constants  t  ,  t  ve  must  discuss  the  structure  of  the  flame 

g*»  ± 

sheet,  where  as  usual  the  coordinate  z  ■  s/ 0  is  appropriate.  The  equations 
for  the  g-perturhat ions  require  their  constancy  through  the  flame  sheet,  so 
that  the  expansions  there  are  written 
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v  ■ T.  *  «v  *  «_1v— *x  • T.  *  8'8.  *  *‘1?i  *•••■ 


(55) 


in  terms  of  vhich 

■  V  -  -  *s.  -  156) 

according  to  the  relations  (3)  and  (4). 

We  first  shov  that  Y  must  he  zero.  Assuming  Y  4  0  and  following 

g*  g* 

the  steps  leading  to  the  structure  problem  (22)  shows  that  a  balanced  equation 
for  can  only  be  obtained  if 


Ta2tn0/(t#  ♦  T  )e 


C57) 


The  corresponding  structure  problem  is 


d^Tj/dz2  ■  Ae 


with  A  *  -  oY  C/km  , 

ge 


(58) 


dtj/dz  -*■  6  ♦  o(l)  as  z  -*•  -«•,  dx^/ds  -»■  o(l)  as  z  -*■  +®.  J 


Clearly  the  boundary  condition  at  s  ■  can  only  be  satisfied  if  A  ■  0. 
Thus  the  assusiption  ^  ^  0  is  incorrect  sad  we  must  have 


Yg#  -  0,  i|#  •  (Y-l)k/2,  Tg#  • 


where  t  ■  (T  *  -1)(y-1)/2. 


(59) 


Balsmcing  the  xj-«quation  with  Yg^  ■  0  leads  to  a  different  condition 


( C/m2kg)exp{eg( t#  -  k£)/T2}  ■  g2/2T^, 


(6o) 


on  g,  namely 
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tbe  T^-structure  being  precisely  that  of  the  combustion  approximation. 
Nov  ve  choose  g  to  be  the  solution  of 


thereby  obtaining  a  positive  guage  function  soley  dependent  on  e  with  the 
property  1  »  g(e)  »0~1.  Then  equation  (60)  becomes 


( C/m2k ) exp {  0g  [  { t#-kt )  /T2  ♦  1]>  ■  g2/2T*, 

&  cL 


from  which  6  disappears  if 


k  -  (t#+T2)/t,  i.e.  ^  -  (t„  +  Tf)g/t. 


This  determines  $  as  a  (linear)  function  of  t». 
o 

We  are  therefore  considering  a  family  of  solutions  for  vhich 

82»2<t..T^)  ,W  <».«*>  r  , 

c - 3? - ’  1-*-  D  - -  *Ip  v%sU) 

mm  * 


.  (64) 


The  family  is  parameterized  by  t#  vhich,  since  is  positive,  cannot 

2 

be  less  than  t#  »  -T  *  so  that  the  present  theory  is  uniformly  valid  in 

ft 

the  interval 


-T  <_  t#  <_  t#0  <  »  for  any  t#e*  (65) 

2 

Elimination  of  t#  and  neglect  of  a  term  0(®g  )  in  the  exponential  yields 


D  ■  ~  —  0  M2  exp[e/(Ta-^t)] 
2^a 


Other  uniformly  valid  families  can  be  obtained  by  taking 


e"0g  ■  ag 


sr 
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vhleh  multiplies  D  by  a  factor  a  .  Because  of  the  exponential,  any 
two  such  families  produce  the  same  D  for  values  of  t#  within  O(0-1g”^) 
of  each  other,  so  corresponding  members  are  identical. 

The  wave  speeds  exhibited  by  {63)  are  vanishingly  small  in  the  limit 
0  -*•  »,  i.e.  much  smaller  than  those  corresponding  to  th»  combustion 
approximation  (which  was  shown  in  Sec.  6  to  make  correct  prediction  for  Mq 
small  but  finite  as  0  -*■  ®. )  Thus  the  solutions  in  this  section  are  very 
slow  deflagrations. 

Now  we  shall  demonstrate  that  the  existence  of  these  other  solutions 
cannot  be  detected  by  the  combustion  approximation,  whose  equations  were 
mentioned  at  the  beginning  of  Sec.  6.  Vith  the  distinguished  limit  and 
the  expansions  of  the  present  section,  it  can  be  shown  that  the  flame-sheet 

perturbation  T  must  again  be  zero  and  that  the  structure  problem  is 

8^ 


(U2)  with 


A  ■  (C/m2kg)exp{0gt„/T^}. 


(68) 


As  usual,  A  is  found  to  be  equal  to  0  /2T  ;  and  that  result  together  with 
(6l)  leads  to  the  analog  of  (62),  i.e. 

.2 


( C/m2k ) exp [ 0g( t  ^/T^+l ) 3  ■  -®r-  , 
a  2T 


(69) 


where  the  choice  (6l)  is  again  made  for  g. 
To  eliminate  0  we  must  set 


•  -  -T2 

'*  ra’ 


(TO) 


thereby  losing  tv  as  a  parameter,  and  then 
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C  -  (B/2Tj)km2,  (71) 

so  that  the  solutions  are  now  parameterized  by  k  ■  >^/g  directly.  The 
corresponding  D  is  (66)  with  t  *  0  which  may  be  rewritten 


a2  2*2 

D  ,  fi.  .a.  A  exp(0/T  ) 
2T  a 


(72) 


2 

on  neglect  of  a  term  O(0g  )  in  the  exponential.  But  this  formula  is  the 

same  as  (44);  we  have  merely  reproduced  the  results  of  Sec.  6,  which  hold 

for  small  M  up  to  0(l)  on  the  0~scale,  in  the  limit  ■  0(g). 

o  o 

Note  that  the  only  difference  between  formulas  (66)  and  (72)  is  t"he  presence 
of  the  term  which  results  from  the  momentum  equation  neglected  in  the 

combustion  approximation. 


8.  Application  of  these  results  for  large  8 

The  mathematical  analysis  provides  an  asymptotic  solution  when  the 
Damkohler  number  has  the  form  (6)  with  C  and  T#  specified.  In  practice, 
however,  D  and  0  (large  but  finite)  are  given  istead  of  C  and  T#, 
which  must  now  be  calculated  from  the  formulas  (47)  and  (48).  Sec.  4  can  then 
be  expected  to  give  an  approximate  solution,  provided  0  is  large.  The 
calculation  of  T#  is  illustrated  in  Fig.  3,  which  graphs  D  as  a  function 
of  T„  in  a  typical  case.  Note  that  D  cannot  exceed  a  certain  value 
corresponding  to  T#  *  (if  a  steady  solution  is  to  exist).  As  6  +  • 

this  upper  bound  increases  indefinitely,  so  that  for  very  large  0  effectively 
any  value  of  D  leads  to  a  steady  deflagration.  As  D  becomes  small, 

T#  ■*  T  antLtKe  solution  obtained  through  the  combustion  approximation  is 


recovered . 


The  procedure  is  particularly  transparent  when  the  heat  release  is 


snail.  Then,  according  to  formula  (36), 


T  -T„ 


*o(T#)  "  (Y.l](T,-l)  * 


(73) 


and  the  restriction  that  0  <  ^  <  requires  that 

—  o  OCJ 


l»CJ 


1  +  B/Y  <  T#  <  T 

* *  a  • 


W) 


Thus 

2  2  T  — T 

D  *  — 6-S— jj-  T—-JT-T-  02exp(0/T#),  (75) 

2(y-1)tI  lT» 

which  clearly  leads  to  a  curve  of  the  form  shown  in  Fig.  3,  the  maximum 
value  of  D  being 

2  2  2 

*  9  u  expted+B/y)].  (76) 

2(1+6/y) 

The  results  in  Sec.  7  show  that  there  is  another  possibility  for 

each  value  of  D.  The  calculation  of  t#  from  the  equation  (6Ub)  for 

large  0  (g  having  been  determined  from  equation  (6l))is  illustrated  in  Fig.  5* 

The  formulas  of  Sec.  6  can  then  be  expected  to  provide  an  approximate  solution. 

Note  that  the  second  solution  exists  whatever  the  value  of  D,  the  wave 

1/2 

speed(on  the  scale  of  g  )  increasing  monotonically  from  0  to  •  as 
D  does  so. 

Thus  we  have  shown  that,  for  a  given  D  in  an  appropriate  range,  there 
are  in  fact  two  distinct  wave  speeds,  one  for  fast  deflagrations  as  predicted 
by  (1*8)  and  illustrated  for  small  heat  release  by  equations  (73)  and  (75) » 
and  one  corresponding  to  very  slow  deflagrations  as  predicted  by  equation  (66). 


In  the  D/n^.M2- plane  both  branches  tend  to  the  single  curve  (UU)  obtained 
by  the  combustion  approximation,  as  M2  0.  A  typical  plot  of  D/m2 
versus  is  shovn  in  Fig.  5. 

The  existence  of  two  wave  speeds  has  been  known  for  some  time, 
albeit  in  association  with  heat  loss  under  the  combustion  approximation. 
Indeed  Spalding  and  Yamlu  (1959)  have  reported  experiments,  involving 
heat  loss,  in  which  two  wave  speeds  were  found.  More  recently,  Buckmaster 
(1976)  described  the  slower  of  the  two  waves  in  the  limit  of  vanishingly 
small  heat  loss.  (The  faster  wave  tends  to  the  familiar  deflagration  wave 
with  T„  ®  T  ,  described  in  Sec.  6  for  0  -*>•». )  However,  a  description  of 
the  very  slow  waves  under  adiabatic  conditions  has  apparently  not  been  given 


before. 
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Appendix  I:  Pr,  Le  ^  1 


The  assumptions  of  Prandtl  and  Lewis  numbers  equal  to  one  were  made  for  sim¬ 
plicity  and  are  not  essential.  Here  we  give  the  changes  that  are  required  when 
we  relax  those  assumptions. 

The  governing  equations  for  general  Pr,  Le  become 

YM2PrdV/ds  -  yM2(V-l)  +  TV-1  -  1  (T7) 

o  o 

dT/ds  +  (Y-l)TV_1dV/ds  -  yd2T/ds2  (78) 

+  Pry(y-1)M2 (dV/ds)2  +  otyYAe"e/T 

dY/ds  =  Le_1d2Y/ds2  -  AYe“6/T  (79) 

Equations  (77)  -  (79)  and  equation  (5)  under  the  assumptions  listed  in  Sec.  2 
are  the  original  equations  that  equations  (1)  -  (4)  were  derived  from.  By  suit¬ 
able  combination  of  equation  ( T7)  and  (78)  one  can  show  that  t  as  defined  by 
(3)  satisfies 

=■  +  (Pr-D(y-l)M2  +  aAYe'e/T  ,  (80) 

ds  ds 

and  from  (80)  and  (79)  H  =  t  +  aY  satisfies 

H  «  1  +  g  +  (y-l)M2/2  +  { (H-x) (Le-1-1)  +  (Pr-1) (y-l)M2V2/2>  .  (8l) 

The  simplification  given  by  the  dual  assumption  that  Pr  *  Le  *  1  is  now  obvi¬ 
ous,  equations  (80)  and  (8l)  reducing  to  (2)  and  (4)  respectively.  Assuming  the 
Prandtl  and  Lewis  numbers  to  be  unity  surely  leads  to  great  simplification  which  j 

motivated  the  form  of  our  original  discussion,  but  generally  the  variables  r  j 

and  II  offer  no  advantage  over  the  original  variables  T  and  Y  which  we  now  1 


will  use  in  the  following. 
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Reactionless  equations  for  the  0(1)  terms  and  0(6  )  perturbations  are 

to  be  subject  to  appropriate  conditions  at  s  ■  -«■>  and  a  set  of  jump  conditions, 
denoted  by 


[u]  -  u+  -  u_  , 


(82) 


derivable  from  a  flame-sheet  analysis.  The  leading  order  conditions  are 


-1. 


[T  ]  =  [V  ]  =  [Y  ]  -  [dV  /ds]  »  [dT  / ds  +  aLe  dY  /ds]  =  0 
0  0  0  0  o  o 


(83) 


T  =  T.  at  s  ■  0  and  Y  ■  0  for  s  >  0 
o  *  o 


The  solution  of  the  leading- order  problem  for  a  given  value  of  Mq  determines  the 


function  T*  =  T*(M  )  as  in  Sec.  4. 


-Is 


Similarily  we  find  conditions  on  the  0(6  )  perturbations  to  be 


-1, 


[Vx]  *  [Tx  +  aLe  ^3=0 

fdVj/ds]  -  [Tjl/yM^PrV,  ,  [d^/ds  +  aLe"1dY1/ds]  -  (1/y-Le) [T^J 


(81) 


The  structure  problem  in  the  flame  sheet  corresponding  to  (22),  now  written  in 


terms  of  T^?  becomes 


d2T 


dz 


1  ~  -  T  /T l 

Y  =  -A  A.  Y^e 


,  T.  +  aLe"1?.  =  (dT  /ds.)z  +  T.  +  aLe^Y. 

1  1  O  +  1+  1+| 


T.  +  (dT  /ds  )z  +  T,  as  z  -  «  ,  T,  +  (dT  /ds,)z  +  T.  as 

1  o  —  1—  -  lo+  1+ 


(8S) 


J 


Z  -v  + 


And  as  before  (85)  can  be  transformed  into  Linan's  problem  (24)  to  be  solved  un¬ 
der  conditions  (26)  with  w  specified  as 
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w  -  -  (dT  / ds , ) / (dT  /ds  -dT  /dsj  ;  (dT  /ds  -dT  /dsj  -  6  .(86) 

Solution  of  Lindn's  problem  then  fixes  the  constants  Yj_,  Y^+  (analogous  to 
t+,  t_  in  Sec.  4), as 


with 


Y_  =  LeT*(n+0_m)/a  ,  Yj  -  Le  T^o^/a 

a  “  l/w[T1+/T^  +  -  £n(82/2ALeTj)] 


(87) 


The  general  problem  as  outlined  is  numerical;  but  the  limit  of  small  heat  re¬ 
lease  yields  explicit  formulas,  as  it  does  in  Sec.  5.  To  all  orders  in  6  the 
solution  for  Yq  is  given  as 

Y  =  Y  (l-eLe  s)  for  s  <  0  ,  Y  -  0  for  s  >  0  .  (88) 

O  O 


Equations  (77)  and  (78).  with  the  reaction  term  set  equal  to  zero  are  to  be  solved 
subject  to  the  required  jump  conditions.  Expanding  Vq  and  Tq  as 


V  -  1  +  Bv*  +...,  T  -  1  +  ST'  +. ... 
o  o 


(89) 


then  v'  and  T'  satisfy 

dv' /ds  *  [ (yM2-1)v'  +  T'J/PryM2  , 
o  o 

d2T'/ds2  ■  [  (y-l)dvVds  ♦  dTVctS  3  /* 


(90) 


subject  to 


T'  0  as  s  -*  -« 


bounded  as 


s  -*•  +  “> 


(91) 


fv'J  =  f  T  *  ]  «  0  ,  [dT'/ds]  -  -1 
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where 


X.,  X  =  {[YM2-l+PrM2]  ±  4YM2-l+PrM2]2  +  4(l-M2)PrYM2}/2PryM2  .(93 

Li  O  O  O  O  O  O  O  ’ 


X^  and  X 2 


have  the  properties: 

X.  ~  1  ,  X0  -  -l/YM2Pr  as  M2  -*•  0  , 

1  2  o  o 

Xx  >  (y-l+Pr) /yPr  ,  X2  --(1-M2)/(Y-1+Pr) 


Setting  T  (0)  equal  to  T*  of  the  definition  (35)  now  shows  that 

(y-l)X2 

t* - r -  *  (95) 

(X2-X1)(1-M^)(yX1-1) 

which  reduces  to  the  formula  (36)  for  Pr  *  1.  As  there  t^  1,  i.e.  -*•  Tg 

as  M2  -+■  0. 
o 

These  formulas  determine  Mq  and  the  combustion  field  to  leading  order  for 

arbitrary  values  of  Pr  and  Le  when  the  heat  release  is  small.  Since  X^ 

2 

is  positive  and  X2  is  negative  for  M*  <  1,  the  result  (92)  shows  that  the 

velocity  v'  increase  monotonically  up  to  the  flame  and  beyond,  ultimately  attain- 

2 

ing  the  value  v'(°°)  *  1/(1-M  ).  On  the  other  hand  the  temperature  increases  up 

2  2 

to  the  flame  but  then  decreases  to  its  final  value  T’(“)  «  (l-YMo)/(l~Mo) .  The 
temperature  and  velocity  profiles  are  therefore  qualitatively  unchanged  from  those 
for  Pr  «*  1.  See  Fig. 2a, b  in  Sec. 4.  However  the  pressure  field,  determined  by  the 
perturbation 
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p'  -  T'  -  v' 


(96) 


is  sensitive  to  the  value  of  the  Prandlt  number.  The  pressure  p*  is  zero  as 

2  2 

s  *  -«®  and  reaches  a  negative  final  value  p'(®)  ■  -yMq/(1-Mo).  In  between,  it 
may  or  may  not  have  a  (positive)  maximum  depending  on  whether  T' (0)  is  greater 
or  less  than  v*(0).  For  Pr  ■  1,  they  are  equal  so  that  the  effect  can  be  ex¬ 
plicitly  exhibited  by  taking  Pr  close  to  one.  Thus  we  evaluate 


p’(0) 


yx2(i-x1) 

U-m^Hx^AjHyXj-i) 


(9T) 


for 


(Pr-1) (y-1)M2  (1-M2) 

XL  -  1 - 2 - 2  and  X2  - 2 — 

vM  +  1  -  M  vM 

o  o  o 


(98) 


to  obtain 


Pf(0)  - 


YM2(Pr-l) 
o _ 

2  2 

YM„+1-M‘ 

O  O 


(99) 


For  Pr  >  1  there  is  a  pressure  spike  at  the  flame,  whereas  for  Pr  <  1  the  de¬ 
crease  is  monotonic  with  s  from  -»  to  +».  Presumably  there  is  a  similar 
effect  when  the  heat  release  is  not  small. 

The  discussion  for  very  slow  deflagrations  found  in  Sec.  7  can  also  be  gen¬ 
eralized  for  Pr,  Le  4  1.  Again  the  starting  point  is  equations  (77)-(T9). 

2 

Considering  the  distinguished  limit  given  by  (50)  and  allowing  Mq  »  kg(0),  the 
analysis  Pr,  Le  t  1  follows  precisely  as  in  Section  7  in  that  we  must  consider 
outer  regions  s  <  0,  s  >  0  and  a  velocity  adjustment  in  the  0(g)  Mach  layer 
behind  the  flame  sheet.  The  detailed  solution  to  0(g)  is  more  complicated  than 
that  presented  In  See. 7f Involving  terms  that  vanish  when  Pr  ■  1,  but  is  obtained 
In  a  straightforward  fashion. 
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Consideration  of  the  flame  structure  requires  as  before  that  1 

thus  equation  (60)  again  holds  with  the  right*hand  side  replaced  with 


82/2T4Le 

a 


Following  the  remaining  arguments  in  Sec.  7  we  find  the  only  change  is 
of  D  parameterized  by  tA  which  becomes 

J-T^N 

0 


8m  V  <'.+V 


2T4Le 

a 


exp  [ 


T.-t*R(0) 


] 


0  and 


(100) 

the  form 


(101) 


It  was  shown  in  Sec.  4  that  for  any  non-zero  Mach  number,  the  combustion  field 


is  comprised  of  a  frozen  region  upstream  and  a  flame  sheet  followed  by  a  frozen 

downstream  region  where  the  temperature  descreases  from  TA.  Thus  reactant  es- 

2 

capes  from  the  flame  in  an  amount  measured  by  the  perturbation  a+  *  ~t+/T*  “ 

aY+/T^  0.  (Note  that  the  partial  burning  exhibited  here  holds  for  all  non- zero 

wave  speeds,  since  Lindn's  problem  controls  the  structure.  In  the  limit 
2 

Mq  -*■  0,  as  in  Sec.  7^Y+  -*  0.)  Because  of  this  non- zero  reactant  fraction  there 
must  be  further  burning  beyond,  albeit  exponentially  weak,  so  as  to  take  the 

reaction  to  completion. 

Consider  the  expansions 


where 


V-V  +  1/0V-  +...,  t-t  +1/0T.+...; 

oo  x  oo  X 

V»  *  V-  ’  -  1  +  6  +  <Y-l>M*/2  ,  Tm  - 


T  -  T  +  1/0T.  +. . . 

V(H4) 

1  +  6  + 


and  the  variable 

-8/T 

00 

n  *  Ae  s  . 


(103) 


"0/T» 

Since  T  <  T.,  and  A  has  the  form  Implied  by  (8),  Ae  °°  is  exponentially 
small  compared  to  0_*.  Hence  the  derivative  in  equation  (1)  is  negligible  to 
all  orders  in  0  ^  and  we  find  from  (1),  (3)  and  (4)  that 

vi  ■  vw*<M>  •  ti  ■  »i  -  fi-WoVt  •  'i  -  -«yi  •  (it*) 


Eliminating  and  In  favor  of  Tj  now  shows  that 


1  *  *  1  »  '  O  m 

~T—  *  -T.expCor.)  with  a  -  —z  ( - 5 - )  , 

dn  1  1  T2  T  -Mv 

00  00  0  ® 

Ti  “  -T*®+  at  n  -  0 

is  the  problem  for  in  n  >  0.  (The  boundary  condition  is  obtained  by  mat¬ 

ching  with  the  burnt  region  immediately  following  the  flame.) 

An  explicit  solution  can  be  given  in  terms  of  the  exponential  integral  function 

E,  (2)  «  f°°e  Cdt/t.  TllttS 
1  z 


A  A  A 

E^(-at^)  -  E^(ao+)  ■  n  for  a  >  0 

Tj  ■  -a+e  '  for  a  ■  0  , 

A  A  A 

E^ari)  -  E^-ao^  -  n  for  a  <  0 

(There  is  no  singularity  of  2  in  (lof)  since  -  M2Vw 


is  positive  for 


(106) 


0  i  Mo  <  Mo<?J  ’ 


as  is  easily  shown.) 
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Captions: 

Fig.  la.  versus  T#  for  Pr  ■  1,  Le  ■  1,y  ■  5/3,  8  ■  .5,  Ta  ■  1.5 

where  .225* 

2 

Fig  lh.  M  versus  t#  for  Pr«l,Le*l,Y*  5/3.  Dashed  line  for 
o  w 

|8»  .09  where  ■  .51.  Solid  line  is  the  asymptotic  result  (34). 

Fig.  2a.  Leading-order  profiles  of  velocity,  temperature,  pressure  and 

mass  fraction  for  Pr  ■  1,  Le  ■  1,y  ■  5/3,8  *  .5,  *  .2  corresponding 

to  T#  ■  1.43. 

Fig.  2b.  Pressure  profiles  for  general  Le,  Y  •  5/3,  *  .5  and  various 

Pr,  showing  a  spike  for  Pr  >  1. 

Fig.  3.  D/m2  versus  T#  for  Pr  ■  1,  La  »  1,  y  *  5/3,  8  *.69,  0  “  50. 

p 

Fig.  4.  D/m  versus  t#  for  Pr  ■  1,  Le  ■  1,  y  ■  5/3,  8  ■  .5,  0  *  50. 

Fig.  5.  M2  versus  D/m2  for  fast  and  very  slow  deflagrations,  Pr  *  1, 

Le  •  1,  y  *  5/3,  3  *  .5,  0  *  50. ^>r  Mack  numbers  af  tht 

upper  it*  4i  of  application,  off  he  core  fusion  «ppK#xtwaf.i»t. 
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The  main  result  is  a  relation  between  the  Mach  number  and  the  flame  tem¬ 
perature,  a  parameter  determined  by  the  Darakohler  number.  One  of  the  features 
of  a  finite  wave  speed  is  the  exponentially  long  tail  in  which  the  small  amount 
of  reactant  escaping  the  flame  front  is  finally  consumed.  Another  is  the  ve¬ 
locity  adjustment  layer  which  forms  behind  the  flame  front  aB  the  Mach  number 
becomes  small,  a  layer  the  combustion  approximation  cannot  detect.  Explicit 
formulas  are  obtained  when  the  heat  release  is  small;  otherwise  numerical  in¬ 
tegration  is  necessary. 

Another  class  of  deflagration  waves  is  uncovered,  characterized  by  Mach 
numbers  that  are  even  smaller  than  those  involved  in  the  combustion  approxi¬ 
mation.  The  existence  of  these  very  slow  deflagrations  has  been  suspected 
before,  but  cannot  be  deduced  from  the  governing  equations  of  that  approxima¬ 
tion. 

For  simplicity  both  the  Prandtl  and  Lewis  numbers  are  taken  to  be  one. 

The  modifications  for  other  values  are  given  in  an  appendix. 


SECURITY  CLASSIFICATION  OF  THIS  FAOIfVton  Data  Kn'aratf) 


